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I.  INTRODUCTION 


New  computational  algorithms  together  with  advances  in  computer  resources 
are  beginning  to  show  promise  as  viable  predictive  capabilities  for  the  aero¬ 
dynamics  of  projectiles.  Recent  papers1’2 3  have  reported  the  development  and 
application  of  the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  computational 
technique  to  predict  the  flow  about  slender  bodies  of  revolution  at  transonic 
speeds.  References  1  and  2  showed  the  technique  to  be  a  viable  computational 
tool  for  predicting  both  external  and  internal  flows  for  spinning  and  non¬ 
spinning  bodies  of  various  geometric  shapes.  These  calculations,  however, 
modeled  the  base  flow  as  an  extended  sting  and  thus  the  base  pressure  and 
reci rculatory  base  flow  were  not  computed.  Only  until  very  recently  with  the 
advent  of  LDV  instrumentation  are  reliable  base  flow  data  becoming  available. 
Therefore,  only  a  very  limited  amount  of  data  are  presently  available  and  in 
general  are  usually  in  terms  of  overall  drag  data  rather  than  specific  base 
pressure  data. 

The  base  flow  problem  at  supersonic  speeds  has  been  an  area  of  extensive 
research.  An  excellent  review  of  base  drag  is  presented  in  Reference  3. 
Recently,  Navier-Stokes  solvers4’5  have  been  used  to  compute  the  aft  end  flow 
field  of  axisymmetric  bodies  at  supersonic  velocities.  At  transonic  speeds,  a 
limited  study  of  the  flow  past  a  boattailed  afterbody  has  been  performed  by 
Chow,  et.  al.6  The  inviscid  transonic  flow  is  solved  by  finite  difference 
calculations  of  the  axisymmetric  potential  equation,  and  the  viscous  flow 
behind  the  base  is  treated  through  integral  formulations.  This  report 


1.  Nietubicz ,  C.J. ,  Pulliam,  T.H.,  and  Stegev,  J.L.,  "Numerical  Solution  of 
the  Azimuthal-Invariant  Thin-Layer  Navier-Stokes  Equations , "  ARBRL-TR- 
02227,  U.S.  Army  Ballistic  Research  Laboratory,  ARRADCOM,  Aberdeen  Proving 
Ground,  MD  21005,  March  1980  (AD  A085716). 

2.  Nietubicz,  C.J.,  "Navier-Stokes  Computations  for  Conventional  and  Hollow 
Projectile  Shapes  at  Transonic  Velocities,"  AIAA  Paper  No .  81-1262,  June 
1981 . 

3.  Sedney,  R.,  "Review  of  Base  Drag,"  U.S.  Army  Ballistic  Research  Labora¬ 
tory,  ARRADCOM,  Report  No.  1337,  Aberdeen  Proving  Ground,  MD  21005, 
October  1966  (AD  808767). 

4.  Diewert,  G.S.,  "A  Computational  Investigation  of  Supersonic  Axisymmetric 
Flow  over  Boattails  Containing  a  Centered  Propulsive  Jet,"  AIAA  Paper  No. 
83-0462,  January  10-13,  1983. 

5.  Weinberg,  B.C.,  McDonald,  H.,  and  Shamroth,  S.J.,  '1 Navier-Stokes  Computa¬ 
tions  of  Aft  End  Flow  Fields, "  Final  Report,  Army  Research  Office  Contract 
DAAG29-79-C-0003,  May  1982. 

6.  Chow,  W.L.,  Bober,  L.J.,  and  Anderson,  B.H.,  "Strong  Interaction  Associat¬ 
ed  with  Transonic  Flow  Past  Boattails,"  AIAA  Jjourngal,  Vol.  13,  No.  1, 

1975,  pp.  112-113. 


7 


describes  a  new  numerical  capability  to  compute  the  flow  field  in  the  base 
region  of  projectiles  at  transonic  speeds  and  thus  the  total  aerodynamic  drag 
can  now  be  computed. 

The  total  drag  for  projectiles  can  be  divided  into  three  components:  (1) 
pressure  (wave)  drag,  (2)  viscous  (skin  friction)  drag,  and  (3)  base  drag. 
For  a  typical  shell  at  M  =  .9,  the  relative  magnitudes  of  the  aerodynamic  drag 
components  are:  (1)  pressure  drag,  20%;  (2)  viscous  drag,  30%;  and  (3)  base 
drag,  50%.  In  order  to  predict  the  total  drag  for  projectiles,  computation  of 

the  full  flow  field  (including  the  base  flow)  must  be  made.  Computation  of 

base  flow  is  especially  important  at  transonic  speeds. 

The  critical  aerodynamic  behavior  of  projectiles  occurs  in  the  transonic 
speed  regime.  This  can  be  attributed  to  the  complex  shock  structure  which 
exists  for  the  projectiles  at  transonic  speeds.  Figure  1  is  a  spark  shadow¬ 
graph  which  shows  the  shock  structure  for  a  typical  projectile  at  M  =  .95, 
a  «  0.  It  also  shows  a  clearly  defined  wake  behind  the  base  of  the  projectile 

devoid  of  any  vortex  shedding.  In  the  present  effort  primary  emphasis  is  on 

the  base  region  flow  field  computations;  however,  the  technique  used  computes 

the  full  flow  field  over  the  projectile  (including  the  base  region).  There¬ 

fore,  all  three  components  of  the  drag  are  computed. 

A  brief  description  of  the  governing  equations  and  the  method  of  solution 
are  given  in  Sections  II  and  III.  A  unique  flow  field  segmentation  procedure 
and  the  implementation  of  boundary  conditions  are  discussed  in  Section  III. 

In  Section  IV  computed  results  are  given  for  transonic  flow  about  a  6-caliber 
secant-ogive-cylinder  shape  for  .9  <  M  <  1.2,  a  =  0.  Velocity  vector  plots 
and  stream  function  contour  plots  are  presented  to  show  the  qualitative 

features  of  the  flow  field  in  the  base  region.  All  three  components  of  drag 
are  obtained.  Base  drag  is  compared  with  experimental  and  semi-empi rical  data 
while  the  total  drag  is  compared  with  available  semi-empirical  data.  The 
encouraging  results  show  that  the  present  computational  technique  can  be  suc¬ 
cessfully  used  to  predict  the  base  region  flow  field  of  projectiles.  The 
results  reported  here  are  for  transonic  speeds,  but  future  computational 
efforts  using  the  same  technique  will  be  directed  at  supersonic  velocities. 


II.  GOVERNING  EQUATIONS 

The  Azimuthal  Invariant  (or  Generalized  Axisymmetric)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  5,  n,  e  can  be  written  as1 

dTq  +  9^E  +  3^G  +  H  =  Re-13^S  (1) 


where 


5  =  5(x,y,z,t)  is  the  longitudinal  coordinate 
n  =  n(y,z,t)  is  the  ci rcumferenti al  coordinate 
C  =  e(x,y,z,t)  is  the  near  normal  coordinate 
t  =  t  is  the  time 
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The  notation  for  the  physical  coordinates  x,  y,  z,  and  the  transformed  coordi- 

A 

nates  K,  n,  ?  are  shown  in  Figure  2.  The  vector  of  dependent  variables  q  and 

AAA 

the  flux  vectors  E,  G,  H  are  given  as 


p 

pU 

, 

pW 

pu 

A  1 

puU+5xP 

puW+?xp 

pv 

-« 

l 

•"O 

II 

UJ 

A 

pvU+5yP 

,  G  =  J  1 

PvW+?yp 

pw 

pwUH  p 

pwW+?zp 

e 

(e+p)U-5tP 

(e+p)W-ctp 

■V  M 

—  — 

L  J 

H  = 


0 

0 

PV[R?(U-5t)  +  Rc(W-ct)] 
-PVR<t>n(V-nt)  -  p/(R*n) 


The  thin  layer  viscous  terms  valid  for  high  Reynolds  number  flow  are  contained 

A  V- 

in  the  vector  S,  where 


S  = 


0 

"('x  +  4  +  'z>uC  +  ("/3)(tx"!  +  Vi  +  ViK 

+  4 *  v\ +  (•j/3)uxuc +  yc +  v^y 
“t«X  +  4  *  SZ)WC  +  (u/3)(Vc  +  Vi  *  Vc)ez 

1UX  +  +  tz)[0.5u(u  +  v  +  w  (j.  +  xPr-1  (t  -  l)_1(a2)s] 

+  (p/3)(cxu  +  Cy v  +  czw)Uxu?  +  eyv?  +  szw^)} 


The  velocities 
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(2) 


U  =  +  5vu  +  E  v  +  5  w 

t  x  y  l 

V  =  nt  +  nxu  +  nyv  +  V 

u  -  «t +  v +  y  +  5zw 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  components  (u,v,w)  are  nondimensional ized  with 
respect  to  a  (the  free  stream  speed  of  sound).  The  density  (p)  is  referenced 
to  and  total  energy  (e)  to  p^a2.  The  local  pressure  is  determined  using 
the  equation  of  state, 

p  =  (y-1 )[e  -  0.5p(u2+v2+w2)]  (3) 


where  y  is  the  ratio  of  specific  heats. 

In  high  Reynolds  number  flows  the  thin-layer  approximation  is  often  used 
because,  due  to  computer  speed  and  storage  limitations,  fine  grid  spacing  can 
only  be  provided  in  one  coordinate  direction.  This  assumption  has  been  vali¬ 
dated  for  flows  about  various  aerodynamic  models  but  still  remains  a  question 
in  the  wake  region  of  bluff  bodies.  The  grid  spacing  available  in  other 
directions  is  usually  too  coarse  to  resolve  the  viscous  terms.  Essentially, 
all  the  viscous  terms  in  the  coordinate  direction  £  and  n  are  neglected  while 
terms  in  the  near  normal  direction  to  the  body  ?  are  retained.  The  thin-layer 
generalized  axisymmetric  equations  (1)  are  obtained  from  the  three  dimensional 
equations  by  making  use  of  two  restrictions:  (i)  all  body  geometries  are  of 
an  axisymmetric  type;  and  (ii)  the  state  variables  and  the  contravariant 
velocities  do  not  vary  in  the  circumferential  direction  (n).  Essentially,  the 
n-derivative  term  in  the  three  dimensional  equations  is  replaced  by  a  source 

term  H  as  it  appears  in  equation  (1).  The  details  can  be  found  in  Reference  1 
and  2. 

Equation  (1)  contains  only  two  spatial  derivatives;  however,  it  retains 
all  three  momentum  equations,  thus  allowing  a  degree  of  generality  over  the 
standard  axisymmetric  equations.  In  particular,  the  circumferential  velocity 
is  not  assumed  to  be  zero,  thus  allowing  computations  for  spinning  projectiles 
or  swirl  flow  to  be  accomplished. 

For  the  computation  of  turbulent  flows  a  turbulence  model  must  be  sup¬ 
plied.  In  the  present  calculations  a  Cebeci-type  two  layer  algebraic  eddy 
viscosity  model  as  modified  by  Baldwin  and  Lomax7  is  used.  In  their  two-layer 
model  the  inner  region  follows  the  Prandtl-Van  Driest  formulation.  Their 
outer  formulation  can  be  used  in  wakes  as  well  as  in  attached  and  separated 


7 .  Baldwin ,  B.S. ,  and  Lomax,  H.,  "Thin-Layer  Approximation  and  Algebraic 
Model  for  Separated  Turbulent  Flows,"  AIAA  Paper  No.  78-257 ,  1978. 
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boundary  layers.  In  both  the  inner  and  outer  formulations  the  distribution  of 
vorticity  is  used  to  determine  length  scales, thereby  avoiding  the  necessity  of 
finding  the  outer  edge  of  the  boundary  layer  (or  wake).  The  magnitude  of  the 
local  vorticity  for  the  axisymmetric  formulation  is  given  by 


O) 


9Wv2  +  /iW  _  iU\2 
3y '  1 3x  3Z ' 


(4) 


III.  NUMERICAL  METHOD 


A.  Computational  Algorithm 

An  implicit  approximate  factorization  finite-difference  scheme  in  delta 
form  is  used  as  described  by  Beam  and  Warming.8 9  An  implicit  method  was  chosen 
because  it  permits  a  time  step  much  greater  than  that  allowed  by  explicit 
schemes.  For  problems  in  which  the  transient  solution  is  not  of  primary 
interest,  this  offers  the  possible  advantage  of  being  able  to  reach  the  steady 
state  solution  faster  than  existing  explicit  schemes. 

The  Beam-Warming  implicit  algorithm  has  been  used  in  various  applica¬ 
tions.1’2’7”10  The  algorithm  can  be  first  or  second  order  accurate  in  time 
and  second  or  fourth  order  accurate  in  space.  The  equations  are  factored 
(spatially  split),  which  reduces  the  solution  process  to  one-dimensional  pro¬ 
blems  at  a  given  time  level.  Central  difference  operators  are  employed  and 
the  algorithm  produces  block  tridiagonal  systems  for  each  space  coordinate. 
The  main  computational  work  is  contained  in  the  solution  of  these  block 
tridiagonal  systems  of  equations. 

B.  Finite  Difference  Equations 

The  resulting  finite  difference  equations,  written  in  delta  form, are 


8.  Beam,  R .,  and  Warming,  R.F.,  "An  Implicit  Factored  Scheme  for  the  Com¬ 
pressible  Navier-Stokes  Equations,"  AIAA  Paper  No.  77-645 ,  June  1977. 

9.  Steger,  J.L.,  "Implicit  Finite  Difference  Simulation  of  Flow  About 
Arbitrary  Geometries  with  Application  to  Airfoils,"  AIAA  Journal,  Vol. 
16,  No.  4,  July  1978,  pp.  679-686. 

10.  Pulliam,  T.H.,  and  Steger,  J.L.,  "On  Implicit  Finite-Difference  Simula¬ 
tions  of  Three-Dimensional  Flow, "  AIAA  Jouvnail,  Vol.  18,  No.  2,  February 
1980,  pp.  159-167. 
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(5) 


(I+h65An-eIJ“1V5A5J)(I+h6cCn-eIJ"1VcA5J 
-  hRe”16?J”1MrJ)  x  (qn+^-qn)  =  -At(<5^En+«5^Gn 
-Re“16?Sn)-AtHn-eEJ“1[(V?A5)2  +  (V^A?)2]Jqn 


Here  h  =  At  because  only  first  order  accuracy  in  the  time  differencing  is 
needed  for  the  steady  state  flows  which  are  considered  here.  This  choice  cor¬ 
responds  to  the  Euler  implicit  time  differencing.  The  6 ’ s  represent  central 
difference  operators,  A  and  V  are  forward  and  backward  difference  operators 

A  A 

^  8E  ^  3G 

respectively.  The  Jacobian  matrices  A  =  — ,  C  =  —  along  with  the  coeffi- 

9q  8q 

A  A 

cient  matrix  M  obtained  from  the  local  time  linearization  of  S  are  described 
in  detail  in  Reference  10.  Fourth  order  explicit  (e^)  and  implicit  (ej)  nu¬ 
merical  dissipation  terms  are  incorporated  into  the  differencing  scheme  to  damp 
high  frequency  growth  and  thus  to  control  the  nonlinear  instabilities.  A 
typical  range  for  the  smoothing  coefficients  is  =  (1  to  5)  At  with 
e  j  —  3eE> 

C.  Flow  Field  Segmentation 


The  objective  of  this  report  is  to  compute  the  full  flow  field  of  a  pro¬ 
jectile  (including  the  base  flow)  at  transonic  speeds.  This  has  been  accom¬ 
plished  by  a  unique  flow  field  segmentation  procedure  which  is  instrumental  in 
the  development  of  the  current  method  for  the  computation  of  base  flow.  This 
segmentation  is  considered  to  be  the  main  contribution  to  the  state  of  art,  as 
far  as  the  computational  algorithm  is  concerned.  The  segmentation  signifi¬ 
cantly  effects  the  treatment  of  boundary  conditions  and  thus  changes  the 
internal  structure  of  the  block  tridiagonal  matrix.  The  details  of  these 
changes  are  described  later. 

Figure  3  is  a  schematic  illustration  of  the  flow  field  segementati on  that 
is  used  to  compute  the  entire  projectile  flow  field  including  the  base  flow. 
It  shows  the  transformation  of  the  physical  domain  into  the  computational 
domain  together  with  the  details  of  the  flow  field  segmentation  procedure  in 
both  domains. 

The  cross  hatched  region  represents  the  projectile.  The  line  BC  is  the 
projectile  base  and  the  region  ABCD  is  the  base  region  or  the  wake.  The  line 
AB  is  a  computational  cut  through  the  physical  wake  region  which  acts  as  a 
repetitive  boundary  in  the  computational  domain.  Implicit  integration  is 
carried  out  in  both  £  and  c  directions.  (See  Figure  3.)  Note  the  presence  of 
the  lines  BC  (the  base)  and  EF  (nose  axis)  in  the  computational  domain.  They 
both  act  as  boundaries  in  the  computational  domain  and  special  care  must  be 
taken  in  inverting  the  block  tridi agonal  matrix  in  the  5  direction. 
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As  a  result  of  the  flow  field  segmentation  procedure  described  earlier 
the  block  tridiagonal  matrix  in  the  K  direction  has  elements  at  J  =  JB,  JB+1 
which  are  treated  as  internal  boundaries  in  the  computational  domain  (J  =  JB 
represents  the  projectile  base  and  J  =  JB+1  is  the  nose  axis).  The  block 
tridiagonal  matrix  in  the  K  direction  takes  the  following  form  (after  setting 
6j  =  0  to  simplify  the  illustration) 


n 


I  a3 

Aq2 

rhs2 

-A2  I  A4 

» 

Aq3 

rhs3 

•  •  • 

•  •  4 

• 

• 

* 

"AjB-2  1  ajb 

• 

• 

0  I  0 

4qJB 

= 

0 

0  I  0 

0 

"aJB+1  1  AJB+3 

•  • 

t 

» 

• 

•  *  • 

• 

m 

“AJMAX-2  1 

• 

AqJMAX-l 

RHSJMAX-l 

— - * 

—  — 

At  A 

Here  A's  denote  the  quantity  a  and  I  is  a  5x5  identity  matrix.  Note  the 

appearance  of  two  uncoupled  block  tridiagonals.  The  rows  at  JB  and  JB+1  are 
particularly  simple  because  boundary  conditions  are  updated  explicitly  at  the 
end  of  inversions.  These  changes  were  easily  implemented  in  a  modular  fashion 
into  an  existing  code  for  projectile  base  flow  computations.  One  simply  fills 
the  block  tridiagonal  matrix  ignoring  the  base  JB  and  the  nose  axis  JB+1. 
Elements  in  these  rows  are  then  overloaded  as  shown  above.  The  flow  field 
segmentation  does  not  affect  the  block  tridiagonal  matrix  in  the  c  direction. 


D.  Implementation  of  Boundary  Conditions 

The  no-slip  boundary  conditions  for  viscous  flow  are  enforced  by  setting 


U  =  V  =  W  =  0 


(7) 


on  the  projectile  surface  except  for  the  base.  At  the  projectile  base  the 
velocity  component  normal  to  the  base  is  set  to  zero,  i.e.  U  =  0,  while  other 
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flow  variables  are  set  equal  to  those  at  a  grid  point  next  to  the  base. 
Future  work  will  be  directed  at  the  implementation  of  viscous  boundary  condi¬ 
tion  at  the  base  to  further  evaluate  this  approximation. 

Along  the  computational  cut  (line  AB)  the  flow  variables  above  and  below 
the  cut  were  simply  averaged  to  determine  the  boundary  conditions  on  the  cut. 
On  the  center  line  of  the  wake  region,  a  symmetry  condition  is  imposed. 


3u 

9z 


=  0 


w  =  0 


(8) 


Free  stream  conditions  are  used  at  the  outer  boundary.  First  order 
extrapolation  for  all  flow  variables  is  used  at  the  downstream  boundary  (lines 
AD  and  AG).  During  transient  calculations,  use  of  a  specified  outflow  pres¬ 
sure  can  give  rise  to  numerical  oscillations  in  the  base  region  flow  field. 
Eventually,  these  grow  and  swamp  the  solution.  This  difficulty  is  avoided  by 
simply  extrapolating  pressure  to  the  downstream  boundary,  which  is  the  proce¬ 
dure  always  used  with  supersonic  outflow.  A  combination  of  extrapolation  and 
symmetry  is  used  at  on  the  nose  axis  (line  EF). 


IV.  RESULTS 

A  series  of  computations  have  been  made  for  the  3 -caliber  (1  caliber  =  1 
max.  body  diameter)  secant-ogive  nose  and  3-caliber  cylinder  shape  shown  in 
Figure  4.  All  the  computations  were  obtained  for  Mach  numbers  .9  <  M  <  1.2 
and  a  =  0.  Limited  base  pressure  measurements  were  made  by  Kayser11  for  this 
projectile  shape  and  are  compared  with  the  computed  results.  The  projectile 
was  supported  by  a  base-mounted  sting  and  meaurements  of  base  pressure  were 
made  at  only  one  location  on  the  base.  These  experiments  were  conducted  at 
Langley  Research  Center  8- foot  Transonic  Pressure  Tunnel.  Computated  base 
pressures  are  also  compared  with  available  semi-empirical12  data.  The  results 
are  presented  in  the  form  of  surface  pressure  distribution,  contour  plots  and 


11.  Kayser ,  L.D. >  "private  communications , "  Ballistic  Research  Laboratory , 
Aberdeen  Proving  Ground MD  21 00 5. 

12.  McCoy,  R.L.,  "McDrag  -  A  Computer  Program  for  Estimating  the  Drag  Coeffi¬ 
cients  of  Projectiles,"  ARBRL-TR-02293,  U.S.  Army  Ballistic  Research 
Laboratory ,  ARRADCOM,  Aberdeen  Proving  Ground ,  MD  21005 ,  February  1981 
(AD  A098110). 
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velocity  vector  plots. 


The  computational  grid  used  for  the  numerical  computations  was  obtained 
from  a  versatile  grid  generator  described  in  Reference  13.  This  program 
allows  arbitrary  grid  point  clustering,  thus  enabling  grid  points  for  the 
projectile  shapes  to  be  clustered  in  the  vicinity  of  the  body  surface.  The 
grid  consists  of  108  points  in  the  longitudinal  direction  and  50  points  in  the 
radial  direction.  The  full  grid  is  shown  in  Figure  5  while  Figure  6  shows  an 
expanded  view  of  the  grid  in  the  vicinity  of  the  projectile.  The  computation¬ 
al  domain  extended  to  4  body  lengths  in  front,  4  body  lengths  in  the  radial 

direction  and  4  body  lengths  benind  the  base  of  the  projectile.  The  grid 

points  in  the  normal  direction  were  exponetially  stretched  away  from  the  sur¬ 
face  with  the  minimum  spacing  at  the  wall  of  .00002D.  This  spacing  locates  at 
least  two  points  within  the  laminar  sublayer. 

The  grid  shown  in  Figure  6  was  generated  in  two  segments.  First,  the 
grid  in  the  outer  region  is  obtained  using  an  elliptic  solver13  for  the  ogive 
portion  and  straight-line  rays  for  the  remaining  portion  which  extends  to  the 
downstream  boundary.  Second,  the  grid  in  the  base  region  is  obtained  simply 
by  extending  the  straight  lines  perpendicul ar  to  line  AB  down  to  the  center 
line  of  symmetry  (line  CD).  An  exponential  stretching  with  the  minimum 
spacing  of  .000020  at  line  AB  is  used.  It  should  be  noted  that  the  same  mini¬ 
mum  spacing  .00002D  is  specified  on  both  sides  of  the  cut, thus  maintaining  a 

smooth  variation  of  grid  across  the  cut.  This  spacing  could,  of  course,  be 

increased  downstream  of  the  base.  The  number  of  grid  points  above  and  below 
line  AB  is  the  same  (50  points)  which  means  that  an  adequate  number  of  points 
is  located  in  the  base  region.  As  can  be  seen  in  Figure  6,  the  grid  points 
are  clustered  near  the  nose-cylinder  junction  and  at  the  projectile  base  where 
appreciable  changes  in  flow  variables  are  expected. 

The  free  stream  Reynolds  number  for  the  series  of  computations  was  fixed 
at  4.5  x  106  based  on  the  total  model  length.  The  computations  are  started 
from  free  stream  conditions  and  marched  in  time  to  obtain  the  steady  state 
solution.  The  initial  calculation  was  made  for  M  =  0.9.  Previous  converged 
solutions  were  then  used  as  starting  conditions  for  additional  Mach  number 
runs  to  achieve  faster  convergence. 

Figures  7  and  8  show  the  distribution  of  the  surface  pressure  coeffi¬ 
cient,  Cp  as  a  function  of  axial  position,  X/D.  Figure  7  shows  the  overall 

view,  whereas  Figure  8  shows  the  distribution  in  the  near  wake  region  of  the 
base.  The  distribution  over  the  projectile  surface  itself  is  shown  in  both 
these  figures.  The  value  of  Cp  beyond  X/D  =  6  is  the  value  of  pressure 

coefficient  along  the  cut  AB.  Both  these  figures  indicate  the  shock  waves 
near  the  nose-cylinder  junction  and  near  the  blunt  base.  Although  not  shown 


13.  Steger,  J.L.,  Nietubics,  C.J.,  and  Heavey ,  K.R. ,  "A  General  Curvilinear 
Grid  Generation  Program  for  Projectile  Configurations ARBRL-MR-03142, 
U.S.  Army  Ballistic  Research  Laboratory ,  ARRADCOM,  Aberdeen  Proving 
Ground,  MD  21005,  October  1981  (AD  A107334) . 
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in  these  figures,  the  pressure  along  the  base  remains  fairly  constant  (within 
±.005  variation  in  Cp  values). 

The  series  of  Figures  9,  10  and  11  show  the  velocity  vector  field  in  the 
base  region  for  M  =  0.9  and  a  =  0.  Each  vector  shows  the  magnitude  and  the 
direction  of  the  velocity  at  that  point.  Figure  9  shows  the  velocity  field  in 
the  entire  base  region.  One  can  see  the  expected  velocity  defect  in  the  far 
wake  region.  Figures  10  and  11  show  the  velocity  field  in  the  vicinity  of  the 
base  (near  wake  region),  the  difference  between  these  two  plots  being  that 
the  former  one  is  stretched  (not  drawn  to  same  scale)  while  the  latter  is 
drawn  with  the  same  scale  in  x  and  y  directions.  Both  figures  clearly  show 
the  reci rcul ati on  region  of  flow  in  the  base  region  and  indicate  a  strong 
shear  layer  as  well. 

The  next  two  Figures  12  and  13  are  stream  function  contour  plots  in  the 
near  wake  region,  again  for  M  =  0.9  and  a  =  0.  Figure  13  is  drawn  to  the  same 
scale  in  x  and  y  while  Figure  12  is  not.  However,  both  of  these  figures  are 
drawn  to  show  the  reci rculation  region  and  the  position  of  the  dividing  stream 
line  as  clearly  as  possible.  They  also  show  the  reattachment  point  which  for 
this  case  is  about  2  calibers  down  from  the  base. 

A  more  critical  check  of  the  computational  results  is  presented  in  Figure 
14  where  the  base  drag  is  plotted  as  a  function  of  Mach  number.  Computational 
results  are  indicated  by  circles,  experimental  results7  by  triangles  and  the 
squares  indicate  the  results  obtained  using  a  semi-empirical  technique  devel¬ 
oped  by  McCoy8.  Base  drag,  as  expected,  increases  as  the  Mach  number 
increases  from  0.9  to  1.2.  The  semi -empi ri cal  technique  shows  generally 
higher  base  drag  when  compared  with  computational  and  experimental  results. 
The  computational  results  predict  the  expected  drag  rise  that  occurs  for  0.9  < 
M  <  1.2.  The  computational  results,  however,  indicate  a  greater  increase  in 
drag  than  predicted  by  either  the  semi -empi ri cal  code  or  the  experimental 
measurements.  The  discrepancy  between  the  numerical  and  the  experimental 
results  should  be  anticipated  since  the  experimental  data  were  obtained  with  a 
sting  attached  to  the  base.  The  sting  has  an  effect  of  weakening  the  recir- 
culatory  flow  in  the  base  region  and  leads  to  higher  base  pressure  and  hence, 
lower  base  drag. 

As  stated  earlier,  the  present  technique  computes  the  entire  flow  field 
about  the  projectile  including  the  base  flow.  Therefore,  in  addition  to 
obtaining  the  base  drag,  the  other  drag  components  (pressure  drag  and  skin 
friction  drag)  are  computed  as  well.  The  three  drag  components  can  then  be 
summed  to  determine  the  total  aerodynamic  drag.  Figure  15  shows  the  variation 
of  pressure  drag  with  Mach  number.  As  Mach  number  is  increased  from  M  =  .9  to 
M  =  1.2,  shock  waves  form  on  the  projectile  surface  and  result  in  increasing 
pressure  drag.  In  Figure  16  skin  friction  drag  has  been  plotted  as  a  function 
of  Mach  number  and  is  shown  to  decrease  as  the  Mach  number  increases  from  .9 
to  1.2.  Figure  17  shows  the  variation  of  total  drag  with  Mach  number  and  is 
compared  with  the  available  semi-empirical  data.  Both  methods  show  the  expec¬ 
ted  drag  rise,  typical  of  transonic  flow  over  projectiles.  The  agreement  is 
shown  to  vary  throughout  the  transonic  regime;  however,  this  represents  the 
first  calculation  of  the  total  aerodynamic  drag  in  a  conceptually  exact  manner 
and  thus  is  very  encouraging. 
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V.  SUMMARY 


A  procedure  has  been  described  in  which  the  Azimuthal-Invariant  (general¬ 
ized  axisymmetric)  thin-layer  Navier-Stokes  code  is  modified  in  such  a  way  as 
to  compute  the  base  flow  field  of  projectiles  at  transonic  speeds. 

The  computed  results  show  the  qualitative  features  of  the  flow  field  in 
the  base  region,  namely  the  recirculation  region,  dividing  stream  line,  reat¬ 
tachment  point,  etc.  Quantitative  comparisons  of  the  base  drag  have  been  made 
with  available  data  for  various  Mach  numbers  in  the  transonic  speed  range. 
These  results  indicate  that  the  present  numerical  technique  can  be  used  to 
calculate  the  base  drag  of  projectiles  at  transonic  speeds. 

The  computed  results  for  this  paper  represent  the  first  application  of 
thin-layer  Navier-Stokes  computational  technique  to  predict  projectile  base 
flow  at  transonic  velocity  using  the  flow  field  segmentation  described  above. 
The  results  indicate  that  this  technique  shows  good  promise  of  providing  a  new 
and  useful  computational  capability  for  exterior  ballistics  of  shells. 

Future  computational  efforts  will  investigate  the  implementation  of 
viscous  boundary  condition  on  the  projectile  base,  improved  grid  resolution, 
alternate  turbulence  models  and  supersonic  flow. 
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Figure  1.  Spark  Shadowgraph  of  Secant-Ogi ve-Cyl inder-Boattail 
Projectile,  M  =  0.95,  a  «  0 
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Figure  2.  Ax i symmetric  Body  and  Coordinate  System 
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Figure  3.  Schematic  Illustration  of  Flow  Field  Segmentation 
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Figure  4.  Model  Geometry 
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Figure  5.  Computational  Grid  for  Flow  Field  Computations 


Figure  6.  Expanded  Grid 
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Figure  7.  Longitudinal  Surface  Pressure  Distribution,  M  =  0.9,  a  -  0 
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Figure  8. 
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Velocity  Vector  Field,  M  =  0.9,  a  =  0 


Figure  10.  Velocity  Vector  Field,  M  =  0.9,  a  = 
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Figure  11.  Velocity  Vector  Field,  M  =  0.9,  a  =  0 
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Figure  13,  Stream  Function  Contours,  M  =  0.9,  a  =  0 
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Figure  14.  Variation  of  Base  Drag  Coefficient  with  Mach  Number,  a 
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Figure  15.  Variation  of  Pressure  Drag  Coefficient  with  Mach  Number,  a 
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Figure  17.  Variation  of  Total  Drag  Coefficient  with  Mach  Number,  a  =  0 
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LIST  OF  SYMBOLS 


a  speed  of  sound 

A  cross  sectional  area  at  the  base 

CD  base  drag  coefficient,  2  D^/p^u^A 

Cp  specific  heat  at  constant  pressure 

Cp  pressure  coefficient,  2(p  -  Po^/p^a* 

D  body  diameter  (57.15mm) 

D b  base  draq 

e  total  energy  per  unit  volume/p^a^ 

A 

q  vector  of  dependant  variables  in  transformed  equations 

A  A 

E,  F  flux  vector  of  transformed  Navier-Stokes  equations 

A 

H  n-invariant  source  vector 

J  Jacobian  of  transformation 

M  Mach  number 

p  pressure/ P^a^ 

Pr  Prandtl  number,  p  C  /< 

oo  n  oo 

R  body  radius 

Re  Reynolds  number,  P^ajl/p^ 

A 

S  viscous  flux  vector 

t  physical  time 

u,v,w  Cartesian  velocity  components/a^ 

U,V,W  Contravariant  velocity  components/a^ 
x,y,z  physical  Cartesian  coordinates 
a  angle  of  attack 

Y  ratio  of  specific  heats 

<  coefficient  of  thermal  conductivity 

p  coefficient  of  viscosity 
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LIST  OF  SYMBOLS  (continued) 

C,n,c  transformed  coordinates  in  axial,  circumferential  and  radial 
di recti ons 

p  densit.y/p^ 

t  transformed  time 

?  circumferential  angle 

Superscri pt 

*  critical  value 

Subscri pt 
b  base 

p  pressure 

v  viscous 

00  free  stream  conditions 
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